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Abstract. It is argued that the consistent description of the transverse- 
momentum spectra, elliptic flow, and the HBT radii in the relativistic heavy-ion 
collisions studied at RHIC may be obtained within the hydrodynamic model if 
one uses the Gaussian profile for the initial energy density in the transverse plane. 
Moreover, we show that the results obtained in the scenario with an early start 
of hydrodynamics (at the proper time to = 0.25 fm) are practically equivalent to 
the results obtained in the model where the hydrodynamics is preceded by the 
free-streaming stage of partons (in the proper time interval 0.25 fm < r < lfm) 
which suddenly equilibrate and with the help of the Landau matching conditions 
are transformed into the hydrodynamic regime (at the proper time tq = 1 fm). 



PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

The problem of the consistent description of various features of the soft hadron 
production in the nucleus-nucleus collisions at the Relativistic Heavy-Ion Collider 
(RHIC) is well known [lj. In particular, the so called RHIC HBT puzzle [TJ d H 0] 
refers to the difficulty of simultaneous description of the hadronic transverse- 
momentum spectra, the elliptic flow coefficient V2, and the Hanbury-Brown-Twiss 
(HBT) interferometry data in various approaches including hydrodynamics [3 [3 

The observable which is particularly difficult to reproduce in the hydrodynamic 
calculations is the ratio of the out to side radius for the pions, -R ut/-Rsidc- In the 
hydrodynamic approach this ratio reaches typically the values around 1.5, while the 
experimental values approach unity (for the average transverse-momentum of the pion 
pair kx ~ 0.5 GeV). The experimental results are usually interpreted as the indication 
for the strong transverse flow and the relatively short emission time. This in turn 
suggests that the equation of state cannot be too soft. In fact, the use of the semi- hard 
equation of state constructed in Ref. |10] , which is based on the smooth interpolation 
between the hadron-gas model and the lattice data, helps to reduce the model values 
of the ratio R ou t/Rside to 1.20-1.25 [11]. At the same time, the transverse-momentum 
spectra and V2 of pions, kaons, and protons are well reproduced. 
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Figure 1. In-plane and out-of-plane sections of the two-dimensional energy- 
density profiles for the centrality class c = 30 — 40% obtained from the Glauber 
model (dashed lines) and the Gaussian parametrization used in Ref. 1141 . Both 
profiles arc normalized to unity, and in addition (x 2 ) = a 2 , (y 2 ) = b 2 . 

Besides the use of the appropriate equation of state, with the temperature 
dependent sound velocity which remains quite substantial at the vicinity of the critical 
temperature T ~ T c , also the use of the complete set of hadronic resonances turns 
out to be responsible for the good description of the data achieved in Ref. 
The inclusion of the feeding from all known hadronic resonances to the yields of 
light particles (such as pions, kaons and protons) is an important feature responsible 
for the success of the thermal models of heavy-ion collisions. It also allows for the 
unified description of the chemical and kinetic freeze-out [12] . The single freeze-out 
assumption shortens the emission time and also helps to reproduce i? ou t /-Rsidc ■ We 
note that in our approach the hadron emission at freeze-out is modeled with the 
Monte Carlo code THERMINATOR [13], which implements the two-particle method of 
the calculation of the correlation function with or without the Coulomb corrections. 

In the recent work [14] . we have shown that the further improvement of the 
description of the data may be achieved if the standard initial conditions, following 
directly from the optical Glauber model, are replaced by the Gaussian energy-density 
profiles of the form 

It is important to emphasize that the Gaussian profiles do originate from the Glauber 
model - they are Gaussian fits to the source distribution determined by the Monte 
Carlo version of the Glauber model [TS] . 

The examples of the in-plane and out-of-plane profiles corresponding to the optical 
Glauber model and the Gaussian parameterization (P) are presented in Fig. [1] We 
note that both types of profiles yield the same values of (x 2 ) = a 2 and (y 2 ) — b 2 . The 
main difference between the two profiles is that the Gaussian profile is steeper in the 
central part. This feature causes the faster development of the transverse flow. Since 
the magnitude of the transverse flow at freeze-out is constrained by the experimental 
data, the faster development of the flow means that the evolution is faster and the 
duration of the emission process is shorter. This is once again the requested feature 
which helps to reproduce correctly the ratio i? ou t/-R s ide - see Fig. O 

Our results presented in Fig. [2] were obtained with the assumption that the 
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Figure 2. The pion HBT radii -R s id c > ^out , fiiongi an d the ratio -Rout/R s ide f° r 
central collisions, shown as the functions of the average momentum of the pair and 
compared to the RHIC Au+Au data. The left panel illustrates our best results 
obtained with the standard Glauber initial conditions while the right panel 
illustrates the results obtained with the Gaussian initial conditions 1 141 . The data 
from | 16| . 



hydrodynamic evolution starts at the proper time tq = 0.25 fm. In the very recent work 
[17] , we have shown that the start of the hydrodynamic evolution may be postponed to 
r = 1 fm if the hydrodynamics is preceded by the free-streaming stage of partons (in 
the proper time interval 0.25 fm < t < 1 fm) which suddenly equilibrate at r = 1 fm. 
In such FS+SE scenario (free-streaming + sudden equilibration, see [HI US]), the 
very fast but delayed equilibration is described with the help of the Landau matching 
conditions for the energy-momentum tensor. Figure [3] illustrates that the results 
obtained in the case when the hydrodynamic evolution starts at tq = 0.25 fm are 
practically indistinguishable from the results obtained within the FS+SE scenario. 
This observation suggests that the consistent description of the soft physics at RHIC 
might be achieved in the approach where the thermalization and the formation of the 
transverse flow happens gradually within 1 fm. 

In conclusion we state that the consistent and uniform description of the soft 
hadronic data at RHIC (including the transverse-momentum spectra, the elliptic flow 
coefficient t>2, the azimuthally sensitive HBT, and non-identical particle correlations 
[TH |2"U1 |2"T] ) may be achieved within the hydrodynamic approach if a proper choice 
of the initial profile and a realistic equation of state are used. Other aspects such as 
the fluctuations of the initial eccentricity [9j[22], inclusion of the resonances, and the 
two-particle method used for the calculation of the correlation function play also an 
important role. Moreover, the start of the hydrodynamic evolution may be delayed 
if the hydrodynamic phase is preceded with the free-streaming stage followed by the 
sudden equilibration. 

We note that a similar conclusion has been also reached recently in Ref. [23] (see 
also [24]) where, however, the fast building of the transverse flow that is required for 
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Figure 3. (Color online) The pion HBT radii -R a id e , Rout , Riong, and the ratio 
Rout/R a ide for central collisions. The darker (lighter) lines describe the results 
with (without) FS+SE. The data from [16]. 



the correct description of the HBT radii is achieved not by the modification of the 
initial conditions but by the very early start of the hydrodynamics (tq = 0.1 fm) and 

the inclusion of the viscous effects. 
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